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Abstract 

The remarkable improvement in the estimates of different cosmological parame- 
ters in recent years has been largely spearheaded by accurate measurements of the 
angular power spectrum of Cosmic Microwave Background (CMB) radiation. This 
has required removal of foreground contamination as well as detector noise bias 
with reliability and precision. Recently, a novel model-independent method for the 
estimation of CMB angular power spectrum from multi-frequency observations has 
been proposed and implemented on the first year WMAP (WMAP-1) data by Saha 
et al. 2006. We review the results from WMAP-1 and also present the new angular 
power spectrum based on three years of the WMAP data (WMAP-3). Previous 
estimates have depended on foreground templates built using extraneous observa- 
tional input to remove foreground contamination. This is the first demonstration 
that the CMB angular spectrum can be reliably estimated with precision from a self 
contained analysis of the WMAP data. The primary product of WMAP are the ob- 
servations of CMB in 10 independent difference assemblies (DA) distributed over 5 
frequency bands that have uncorrelated noise. Our method utilizes maximum infor- 
mation available within WMAP data by linearly combining DA maps from different 
frequencies to remove foregrounds and estimating the power spectrum from the 24 
cross power spectra of clean maps that have independent noise. An important merit 
of the method is that the expected residual power from unresolved point sources is 
significantly tempered to a constant offset at large multipoles (in contrast to the 
~ I 2 contribution expected from a Poisson distribution) leading to a small correc- 
tion at large multipoles. Hence, the power spectrum estimates are less susceptible 
to uncertainties in the model of point sources. 
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1 Introduction 



Remarkable progress in cosmology has been made due to the measurements 
of the anisotropy in the cosmic microwave background (CMB) over the past 
decade. The extraction of the angular power spectrum of the CMB anisotropy 
is complicated by foreground emission within our galaxy and extragalactic 
radio sources, as well, as the detector noise |l[ 0). It is established that the 
CMB follows a blackbody distribution to high accuracy (0). Hence, foreground 
emissions may be removed by exploiting the fact that their contributions in 
different spectral bands are considerably different while the CMB power spec- 
trum is same in all the bands (0; H; S 0). Different approaches to foreground 
removal have been proposed in the literature (|l|; L=T m |9j llOr ) . 

The Wilkinson Microwave Anisotropy Probe (WMAP) observes in 5 frequency 
bands at 23 GHz (K), 33 GHz (Ka), 41 GHz (Q), 61 GHz (V) and 94 GHz (W). 
In the first data release, the WMAP team removed the galactic foreground 
signal using a template fitting method based on a model of synchrotron, free 
free and dust emission in our galaxy (J6J). The sky map around the galactic 
plane and around known extragalactic point sources were masked out and the 
CMB power spectrum was then obtained from cross power spectra of inde- 
pendent difference assemblies in the 41 GHz, 61 GHz and 94 GHz foreground 
cleaned maps fl2l). 

A model independent removal of foregrounds has been proposed in the lit- 
erature (0). The method has also been implemented on the WMAP data in 
order to create a foreground cleaned map (flih. The main advantage of this 
method is that it does not make any additional assumptions regarding the na- 
ture of the foregrounds. Furthermore, the procedure is computationally fast. 
The foreground emissions are removed by combining the five different WMAP 
bands by weights which depend both on the angular scale and on the location 
in the sky (divided into regions based on 'cleanliness'). However, the power 
spectrum recovered from the single foreground cleaned map has severe excess 
power at large multipole moments due to amplification of detector noise bias 
in the autocorrelation power spectrum beyond the beam resolution. 

The prime objective of our method 0) is to adapt and extend the foreground 
cleaning approach to power spectrum estimation by retaining the ability to 
remove detector noise bias exploiting the fact that it is uncorrelated among the 
different Difference Assemblies (DA) {3 El). The WMAP data uses 10 DA's 
f 15l ; 16l : 17; HI), one each for K and Ka bands, two for Q band, two for V band 



and four for W band. We label these as K, Ka, Ql, Q2,V1, V2, Wl, W2, W3, 
and W4 respectively. We eliminate the detector noise bias using cross power 
spectra and provide a model independent extraction of CMB power spectrum 
from WMAP first year data. So far, only the three highest frequency channels 
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observed by WMAP have been used to extract CMB power spectrum and the 
foreground removal has used foreground templates based on extrapolated flux 
from measurements at frequencies far removed from observational frequencies 
of WMAP (Il2t Il9t I20T ). We present a more general procedure where we use 
observations from all the five frequency channels of WMAP and do not use 
any observational input extraneous to the WMAP data set. 



2 Methodology 



The entire method can be neatly split into three distinct sections - 

• foreground cleaning, 

• power spectrum estimation, and, 

• correcting for known systematics. 

We discuss in some detail the systematic correction for residual power due to 
unresolved sources 1 . A merit of our approach is that the expected residual 
power in the cross spectra from combined maps has a flat dependence with 
increasing multipoles - a strongly tempered dependence given that the con- 
tribution in individual DA maps scales as I 2 . (This feature holds promise for 
experiments at higher angular resolution.) 

The three parts of the current methodology are logically modular and it is 
possible to modify and improve each one relatively independent of the other. 
Specifically, we have employed the Pseudo-C; (MASTER (|23h) power spectrum 
estimation from the foreground cleaned maps which is fast but sub-optimal for 



the low I (in practice). More appropriate approaches have been proposed ((24; 
25l : 26h . We are exploring improvements on this front. 



2.1 Foreground Cleaning 



The foreground cleaning stage of our method is adapted from Tegmark & 
Efstathiou (1996) and Tegmark et al. (2003). In Tegmark et al. (2003), a fore- 
ground cleaned map is obtained by linearly combining 5 maps corresponding 
to one each for the different WMAP frequency channels. For the Q, V and W 
frequency channels, where more than one maps were available, an averaged 
map was used. However, averaging over the DA maps in a given frequency 
channel precludes any possibility of removing detector noise bias using cross 
correlation. In our method we linearly combine maps corresponding to a set 
of 4 DA maps at different frequencies. We treat K and Ka maps effectively 

1 Study of other effects, such as beam non-circularity is in progress |2ll : 
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(K,KA)+Q1+V1+W12= 
(K,KA)+Q1+V1+W13= 
(K,KA)+Q1+V1+W14= 
(K,KA)+Q1+V1+W23= 
(K,KA)+Q1+V1+W24= 
(K,KA)+Q1+V1+W34= 
(K,KA)+Q2+V2+W12= 
(K,KA)+Q2+V2+W13= 
(K,KA)+Q2+V2+W14= 
(K,KA)+Q2+V2+W23= 
(K,KA)+Q2+V2+W24= 
(K,KA)+Q2+V2+W34= 



(Cl.CAl) 

(C2,CA2) 

(C3.CA3) 

(C4.CA4) 

(C5.CA5) 

(C6,CA6) 

(C7,CA7) 

(C8,CA8) 

(C9,CA9) 

(C10,CA10) 

(Cll.CAll) 

(C12.CA12) 



(K,KA)+Q1+V2+W12= 
(K,KA)+Q1+V2+W13= 
(K,KA)+Q1+V2+W14= 
(K,KA)+Q1+V2+W23= 
(K,KA)+Q1+V2+W24= 
(K,KA)+Q1+V2+W34= 
(K,KA)+Q2+V1+W12= 
(K,KA)+Q2+V1+W13= 
(K,KA)+Q2+V1+W14= 
(K,KA)+Q2+V1+W23= 
(K,KA)+Q2+V1+W24= 
(K,KA)+Q2+V1+W34= 



(C13,CA13) 
(C14.0A14) 
(C15,CA15) 
(C16.CA16) 
(C17.CA17) 
(C18.CA18) 
(C19.CA19) 
(C20,CA20) 
(C21.CA21) 
(C22.CA22) 
(C23.CA23) 
(C24,CA24) 



Table 1 

List of the 48 possible combinations of the DA maps that lead to 48 cleaned maps 
(Ci, CAi), i=l,2,..,24.. 

as the observation of CMB in two different DA (at the lowest frequencies). 
Therefore we use K and Ka maps in separate combinations. In case of W 
band 4 DA maps are available. We simply form an averaged map taking two 
of them at a time and form effectively 6 DA maps corresponding to W band. 
Wij represents simply an averaged map obtained from the i th and j th DA of 
W band 2 . In table 1 we list all the 48 possible linear combinations of the 
DA maps that lead to 'cleaned' maps, Ci and CAi's, where i = 1, 2, . . . , 24. 

The goal is to determine a set of optimal weights for each multipole, 



[W, 



(wl,wf,wf,wf 



that specifies the linear combination of 4 DA in the combination, which lead 
to a 'cleaned' map with minimal power 



i=4 i 



afJr = E^' (2) 



=i D l 



where a\ m is spherical harmonic transform of map and B\ is the appropriate 
beam function for the channel i. The condition that the achromatic CMB 
signal remains untouched during cleaning is encoded as the constraint 

[W I ][e] = [e] T [W] T = l, (3) 

2 Other variations are possible and some explored. We defer a discussion to a more 
detailed publication (27). 
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where [e] = (1, 1, 1, 1) T is a 4 x 1 column vector with unit elements. The 
constraint implies that the total power in the cleaned map 

c Clean = [Q] ^jT = ^5 + [C jR)] ^jT (4) 

where Cf is the CMB power spectrum and the last term is the unwanted (non- 
cosmic, frequency dependent) contribution from the foreground and other con- 
taminants. 

The optimum weights to combine 4 different frequency channels such as to 
minimize the second term on the rhs of eq. 4 subject to the constraint that 
CMB is untouched (eq. 3) is readily obtained as ((2; EH) 

[W I ] = [e] T [C I ]- 1 /([e] T [C I ]- 1 [e]). (5) 



Here the matrix 



1 m=l n* n j * 

[CJ EE C? = — — Y =^ . (6) 

[ IJ 1 21 + 1 ^ B\B\ 1 ; 



In practice, we smooth all the elements of the Ci using a moving average 
window over A/ = 11 before deconvolving by the beam function. This avoids 
the possibility of an occasional singular C; that cannot be inverted. The entire 
cleaning procedure is automated and takes approximately 3 hours on a 16 
alpha processor machine to get the 48 cleaned maps. In all the 48 maps some 
residual foreground contamination is visibly present along a small narrow strip 
on the galactic plane. For the angular power estimation that follows, the Kp2 
mask employed suffices to mask the contaminated region in all the 48 maps. 
In ongoing work, we assess the quality of fore grou nd cleaning in these maps 



using the Bipolar power spectrum method (|28t l29l : l30t l3ll ). 

2.2 Power Spectrum Estimation 



Independent estimates of the power spectrum can be obtained by cross corre- 
lation power spectra 

cf = {2i + i)- 1 Y.o i iM m y (7) 



of pairs of cleaned maps Ci & Cj chosen such that they share no common 
DA. This is ensures a zero noise bias of the estimates. A Kp2 mask is applied 
to the maps to mask out foreground contaminated regions and accounted 
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Fig. 1. The final angular power spectrum from WMAP-1 (red) and WMAP-3 (blue 
& green) data are shown in the upper panel. The lower panel shows the correction 
made for the residual power from unresolved point sources based on the point source 
models published by WMAP team for the corresponding year. The green curve in 
the panels correspond to the WMAP-3 spectrum (upper), and the corresponding 
correction (lower) obtained using a revised point source model suggested by Huf- 
fenberger et al. 2006. Note the tempered behavior of the point source correction in 
our method (also see fig. 4). The multipole range is log scale for I < 100, and linear, 
thereafter. 



for by de-biasing the pseudo-C/ estimate using the coupling (bias) matrix 
corresponding to the Kp2 mask, appropriate circularized beam transform and 
pixel window (jish . The left panel of figure 2 plots and lists the 24 cross power 
spectra for which the noise bias is zero and closely match each other for I < 
540. (This is in contrast to pre-point source corrected cross power spectra in 
WMAP-1 power spectrum estimation shown in the right panel - a merit of 
our method that we discuss in the next section.) 

An 'Uniform average' power spectrum is then obtained by combining the 24 
cross power spectra with equal weights 3 . The power spectrum is then cor- 
rected for residual power from unresolved point sources as described in the 
next subsection. The error bars have been estimated with comprehensive sim- 
ulations described in §2.4. The final power spectrum is binned in the same 
manner as the WMAP's published result for ease of comparison. 



3 There exists the additional freedom to choose optimal weights for combining the 
24 cross-power spectra. 
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24 Cross-Power Spectra 
C1 ®CA12 C13 0CA24 
C2®CA11 C14 0CA23 
C3®CA10 C15 0CA22 
C4®CA9 C16 0CA21 
C5 ® CA8 
C6 ® CA7 
C7 ® CA6 
C8 ® CA5 
C9 ® CA4 



C10®CA3 C22 0CA15 
C11 ®CA2 C230CA14 
C12®CA1 C24 0CA13 




Fig. 2. Left: Plot of the 24 individual cross power spectra corresponding to the 
cross correlations listed in this figure are show very small dispersion for I < 540. 
The average power spectrum is plotted in red line and blue error bars. Right: Plot 
of the 28 cross-power spectra in the WMAP-1 team analysis prior to point source 
corrections show a large scatter at large /. The multipole range is log scale for 
I < 100, and linear, thereafter. 

The power spectra from our independent analyzes is entirely consistent with 
the WMAP-1 (HI 13) and WMAP-3 (0 0) team results. We find a sup- 



pression of power in the quadrupole and octopole moments consistent with 
WMAP-1 published result. However, in WMAP-1 our quadrupole moment 
(U6fiK 2 ) is a little larger than WMAP-1 result (123/iiT 2 ) and Octopole (455/xK 2 ) 
is less than WMAP-1 result (611fiK 2 ). Our WMAP-1 spectrum does not show 
the 'bite' like feature present in WMAP's power spectrum at the first acoustic 
peak reported by WMAP-1 (jlil b In the three year results, the differences at 
low / < 10 between WMAP result and our is primarily due to the fact that 
WMAP report the maximum likelihood values for I < 10 (and we employ 
pseudo-Q estimation). At large /, our results are consistent with other esti- 
mates as shown in a comprehensive reanalysis of WMAP-3 carried out jointly 
with five independent, international CMB analysis groups (see MASTERint 
in Ref. (0)). 

A quadratic fit to the peaks and troughs of the binned WMAP-1 [WMAP-3] 
power spectrum shown in the figure 3, locates the first acoustic peak at / = 
219.8±0.8 [219.9 ± 0.8] with amplitude AT) = 74.1 ±0.3 [744 ± 0.3]fiK , 
the second acoustic peak at I = 544±17 [539.5 ± 3.8} with amplitude AT; = 
48.3±1.2 [49.4 ±0.4}fxK and the first trough at Z = 419.2±5.6 [417.7 ± 3.2] 
with amplitude AT; = 41.7 ± 1 [41 .3 ± 0.6][lK. 



2. 3 Residual unresolved point source correction 



The residual power contamination in the 'Uniform average' power spectrum 
from unresolved point sources can be estimated by running through our anal- 
ysis the frequency dependent power obtained from the same source model 
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Fig. 3. The red line shows the quadratic function fitted to the peaks and troughs 
of the final binned angular power spectrum for WMAP-1 data (left) and WMAP-3 
(right). SI, S2 and S3 are the 3 different sections of the spectrum on which fits are 
performed individually. The power spectrum estimates with error bars and multipole 
band are represented as the shaded boxes. 

used by WMAP team to correct for this contaminant (Jlii WMAP has a 
fitted a simple model with a constant spectral index over the entire sky for 
the residual power in cross-spectra c\ P of frequency channel % and j from 
unresolved point sources as (fj) 



(8) 



The amplitude A and slope (3 ~ —2 has been obtained by extrapolation from 
the fluxes and spectra of 208 resolved point sources in the WMAP-1 (J6) and 
300 point sources in WMAP-3 Q. 

The contribution from each region a = 1,2,. ..,9 to the expected residual 
from unresolved point sources in our cross-power spectrum of cleaned maps 
i and j which are each linear combinations of 4 DA's at different frequencies 
(eq. 2) is given by 



where w\, (i) denotes the weight vector for the cleaned map i and v$ is a fiducial 
frequency. The index pair (ij) are effectively a single index limited to the 24 
'allowed' cross-spectra. The total contribution is added up as 

24 

Cf S > = (1/24) E E E MtfctfSW (10) 

(ij)=l a I' 



where the coupling matrix for the individual regions where a = 1, 2, 9. 
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Fig. 4. Comparison of the residual power in WMAP due to unresolved point source 
contamination (black) in the final averaged cross power spectra in our method with 
the actual level point source contamination in each the 5 frequency channels. Note 
that the level of residual contamination is in between the cleanest (highest) two 
frequency channels (V & W). More importantly, the residual point source correction 
is constant with multipoles whereas the actual level of contamination increases oc I 2 
with multipoles. This implies that our method is efficient in canceling out the point 
source contamination at high resolution. 

We use the above expression to estimate the unresolved residual point source 
contamination in the final cross power spectrum. We note in passing that this 
estimation is entirely based upon the assumption that the residual unresolved 
point source contamination is statistically isotropic over the sky. The residual 
power from unresolved point sources in WMAP-1 [WMAP-3] is a constant 
offset of ~ 140 [100}fiK 2 for I > 400 (and negligible at lower I). As seem 
in figure 4, this residual is much less than actual point source contamination 
in Q, KA or K band and intermediate between V and W band point source 
contamination. It is noteworthy that the method significantly tempers the 
point source residual at large I that otherwise is oc I 2 in each map. 



2.4 Error Estimate on the Power Spectrum 



The errors on the final power spectrum are computed from Monte Carlo sim- 
ulations of CMB maps for every DA each with a realization of the WMAP 
noise and common signal and diffuse foreground contamination. For WMAP- 
1, we use the 110 random noise maps made available at the LAMBDA data 
archive. For WMAP-3 year analysis random noise map corresponding to each 
DA are generated by first sampling a Gaussian distribution with unit variance. 
In the final step we multiply each Gaussian variable by the number ao/)J~Np 
to form realistic noise maps. Here o~o is the noise per observation for the DA 
under consideration and \Jn~ p is the effective number of observations at each 
pixel. We use the publicly available Planck Sky Model to simulate the con- 
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tamination from the diffuse galactic (synchrotron, thermal dust and free-free) 
emission at the WMAP frequencies. The CMB maps were smoothed by the 
beam function appropriate for each WMAP's detector. The set of realistic DA 
maps with noise, foreground and CMB are then passed through the cleaning 
pipeline. The common CMB signal in all the maps was based on a realization 
of the WMAP 'power law' best fit A-CDM model (|34T). Averaging over the 
power spectra from the simulations we recover the model power spectrum, 
but for a hint of bias towards lower values in the low I moments. For I = 2 and 
/ = 3 the bias is —27.4% and —13.8% respectively. However this bias become 
negligible at higher I, e.g. at I = 22, it is only —0.8%. 

The standard deviation obtained from the diagonal elements of the covariance 
matrix is used as the error bars on the C/'s obtained from the data 4 . The 
covariance matrix of the binned power spectrum is largely diagonal. 



3 Conclusion 

The rapid improvement in the sensitivity and resolution of the CMB experi- 
ments has posed increasingly stringent requirements on the level of separation 
and removal of the foreground contaminants. Standard approaches to fore- 
ground removal, usually attempt to incorporate the extra information about 
the amplitude, spatial structure and distribution of foregrounds available at 
other frequencies, in constructing a foreground template at the frequencies of 
the CMB measurements. These approaches could be susceptible to uncertain- 
ties and inadequacies of modeling involved in extrapolating from the frequency 
of observation to CMB observations. 

We carry out an estimation of the CMB power spectrum from the WMAP data 
that is independent of foreground model and evades these uncertainties. The 
novelty is to make clean maps from the difference assemblies and exploit the 
lack of noise correlation between the independent channels to eliminate noise 
bias. This is the first demonstration that the angular power spectrum of CMB 
anisotropy can be reliably estimated with precision solely from the WMAP data 
(difference assembly maps) without recourse to any external data. Our work 
is a clear demonstration that the blind approach to foreground cleaning is 
comparable in efficiency to that from template fitting methods and certainly 
adequate for a reliable estimation of the angular power spectrum. Moreover, 
this method enjoys the advantage of eliminating the uncertainties of modeling 
required to create the foreground templates that estimate the contamination at 
the CMB dominated frequencies by extrapolating from the observed emission 

4 The beam uncertainty is not included here, but is deferred to future work where 
we also incorporate non-circular beam corrections |2lh. 
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flux at very different frequencies where the foregrounds dominate. 

Further, the tempered, flat, large I behavior of residual from unresolved point 
sources in this method holds promise for experiment with finer angular reso- 
lution. The understanding of polarized foreground contamination in CMB po- 
larization maps is rather scarce. Hence modeling uncertainties could dominate 
the systematics error budget of conventional foreground cleaning. The blind 
approach extended to estimating polarization spectra after cleaning CMB po- 
larization maps could prove to be particularly advantageous. 
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